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Abstract 



We consider arbitrary splits of field operators into two parts; ij) = tp + +^~ , 
and use the corresponding definition of normal ordering introduced in [1]. In 
this case the normal ordered products and contractions have none of the 
special symmetry properties assumed in existing proofs of Wick's theorem. 
Despite this, we prove that Wick's theorem still holds in its usual form as 
long as the contraction is a c-number. Wick's theorem is thus shown to be 
much more general than existing derivations suggest, and we discuss possible 
simplifying applications of this result. 
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I. INTRODUCTION 



Wick's theorem is at the core of perturbative calculations in the canonical operator 
approach to quantum field theory: on expansion of the ^-matrix, it is used to relate time 
ordered products of field operators to normal ordered products and contractions [2] . Within 
the traditional definition of normal ordered products, the interaction picture field operators, 
ip(x), are split into positive ^ 0+S) and negative -?// 0- * 1 energy waves containing annihilation 
and creation operators respectively; 

^(x)=^ 0+ \x)+^°-\x). (1.1) 

Normal ordering, 7V(°)[«], then places the annihilation (creation) operators to the right (left), 
otherwise leaving the order of the operators unchanged. Hence the vacuum expectation 
values of normal ordered products are guaranteed to vanish, 

<0|A^<°>[.]|0> = 0. (1.2) 

Calculations therefore proceed in a straightforward manner by simply taking the vacuum 
expectation value of the operator identity called Wick's theorem. 

In the literature, Wick's theorem has been proved either through the use of a generating 
functional [3-5], or by the manipulation of products of operators [2,6-8]. In both cases, 
an essential ingredient of the proof is the fact that the normal ordered products defined 
above are symmetric in the following sense. Let us denote the normal ordered product of n 
(possibly different) field operators := ip(xi) evaluated at different space time points Xi by 



Then since 



C« :=Jv(0) [^-il' (1-3) 

^(°+)(a;),^( 0+ )(a;')l = \ip {0 ~Hx), -0 (o_) «)l = Va ^' ( where ° = + 1 for 
bosonic operators, and a — — 1 for fermionic ones), it follows that 

= (-l) P <t.,« Vn > 0. (1.4) 

Here a,b,...,t is any permutation of l,2,...,n and p is the number of times pairs of 
fermionic fields are interchanged in moving from an order 1,2, ... ,n to an a,b, . . . ,t or- 
der of fields. Normal ordered products satisfying (1.4) are said to be symmetric. The reason 
why this symmetry is central to the existing derivations of Wick's theorem will be explained 
in detail in section 3 below where we work with the generating functional approach. How- 
ever, the following comments may help illustrate the reason. Recall that Wick's theorem 
relates time ordered products, normal ordered products and contractions. One way to prove 
this theorem [2,7] is to consider initially a given order of times. The resulting (single) con- 
tribution from the time ordered product may then shown to be equivalent to a combination 
of normal ordered products and contractions. If these normal ordered products and con- 
tractions are symmetric, then each contribution to the time ordered product takes the same 
form, and thus Wick's theorem is proved [2,7]. However, if they are not symmetric this is 
no longer the case (see section 3). 

In general, however, normal ordering with respect to annihilation and creation operators 
does not simplify calculations. As an example, consider how, in the presence of symmetry 
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breaking, one must normal order with respect to new annihilation and creation operators 
related by a Bogoliubov transformation to the ones appearing in the field [6,7]. This is 
equivalent to choosing to split the field in a different way from (1.1). A more recent example 
[1] was provided by a system in thermal equilibrium at a temperature T. There the thermal 
expectation value of a normal ordered product of creation and annihilation operators (af 
and a respectively) no longer vanishes in the way that it did at T = 0: 

<iV (0) [a,a t ] > = n be ,n /d ^0. (1.5) 

(Here n& e , nfd are the Bose-Einstein and Fermi-Dirac distributions respectively.) As a result, 
perturbative calculations at finite T would not seem to be simplified by the use of Wick's 
theorem. This point has caused some confusion in the literature [7,9,10] and it was addressed 
in [1]. There it was shown that Wick's theorem could be used in the usual way at finite 
temperature provided one split the field in a different way from the usual one into creation 
and annihilation operators (1.1). 

Thus in general one must consider arbitrary splits of the field operators. The precise 
split used will reflect the physical environment of the problem if it is to simplify practical 
calculations in the usual way. We denote the arbitrary split by 

tj}{x) = ip + {x) +i/>~(x), (1.6) 

where ip + (x) and ip~(x) are unspecified in this paper, so that we can accommodate any 
physical problem. We still call the two parts 'positive' and 'negative' even though the split is 
general and no longer necessarily into positive and negative energy waves. The corresponding 
(natural) generalisation of normal ordering is then expressed in terms of ip + (x) and ip~(x), 
and is defined to be 

• Normal ordering of products of fields, iV[«], places all the 'positive' parts to the right 
and the 'negative' parts to the left, otherwise leaving the order of the fields unchanged. 
A sign change occurs whenever two fermion field operators are interchanged. 

Note that this reduces to the usual definition, iV^[«], when the split is into the annihilation 
and creation operators appearing in the field. We have used the (0) superscript to denote 
this traditional and very special case. So, for example, in the two-point case one has 

N fo^ 2 ] := N h2 = + + ^2^t + i>ii>2- (1-7) 

Such a generalisation of normal ordering was successfully used for thermal field theory in 
[1], and is equivalent to discussions of normal ordering involving Bogoliubov transformations 
and symmetry breaking [6,7]. Contractions are then defined in the usual way through the 
two-point time ordered products which, being independent of the split, take the canonical 
form; 

Ti, 2 := T [^2] = ^ 2 0(fi - h) + <7^i0(* 2 - fi). (T2) (1.8) 

(In general the 9 functions may be defined for complex times and hence be contour dependent 
- as for example in Euclidean quantum field theory and thermal field theory [3,11]. Thus 
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the time ordered products are also contour dependent. However, the analysis presented here 
holds regardless of contour.) Using (1.7) and (1.8) the contraction D 12 is: 



D 1>2 = D[i/; 1 ,i/> 2 ] :=T 1>2 -N 1>2 

= d(h-t 2 ) 

L J (7 

- o(t 2 - h) { Uf, ^ 2 + l + k~, ^t] + k~ ^2~1 } • 



(1.9) 



1.10) 



Now the key point is that for such arbitrary splits, the contraction and two-point normal 
product are not generally symmetric: from (1.7) and (1.10) the condition for symmetry is 
that for all times ti,t 2 



D 12 = aD 2A , N lj2 = aN 2A 



+ 



^1 ,^2 



o, 



(1.11) 



which is not necessarily satisfied. If these two-point relations are true, then all higher order 
normal ordered products are symmetric in the sense of (1.4). This is only easy to prove 
if one has already proved that Wick's theorem holds independent of the symmetry of all 
normal ordered products! 

It is conceivable that in some situations, for example systems out of equilibrium, calcu- 
lations may be simplified by the use of a split for which the normal ordered products are 
not symmetric (i.e. (1.11) is not satisfied). Indeed, out of equilibrium it may be possible to 
find a split such that (1.11) is true for specific times, but without time-translation symmetry 
it is doubtful that this split will ensure (1.11) will be true for all times. Wick's theorem, 
which relates normal ordered products (1.7), time ordered products (1.8), and contractions 
(1.10), has only been proved for the special case where the normal ordered products and 
contractions are symmetric. The question this paper asks then is: Does the usual form of 
Wick's theorem still hold when normal ordered products are not symmetric? 

Our answer is yes. Wick's theorem is in fact more general than the existing proofs 
suggest, namely: if the contraction is a c-number (and hence the split is linear in creation 
and annihilation operators), Wick's theorem holds for both symmetric and non-symmetric 
normal ordered products. Possible practical applications of this result, in particular to time 
dependent non-equilibrium systems, are discussed in the conclusions. 

In section 2 we derive the generating functional for Wick's theorem in terms of the 
generalised definition of normal ordering. We stress the places in which the non-symmetry 
of products means that our derivation differs from others existing in the literature. In section 
3 we use the generating functional to obtain Wick's theorem. For symmetric products the 
proof is immediate - for non-symmetric ones it is considerably more complicated. We finally 
summarise our work and also comment on a claim made in the original paper on thermal 
field theory by Matsubara [9]. Applications for our work are then discussed. 

II. THE GENERATING FUNCTIONAL FOR WICK'S THEOREM 



In this section we derive the generating functional for Wick's Theorem - that is, we prove 
the identity 
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T 



x 



exp{— i J d 4 x j(x)ip(x)} — N exp{— % J d A x j(x)ip(x)} 
exp{-^ J d 4 x d 4 y j(x)D[ip(x)ip(y)]j(y)}. 



(2.1) 



Here / d 4 x = J dt J d 3 x, and the precise range of integration does not affect our arguments. 
For simplicity, we just consider bosonic fields. (Comments regarding fermionic fields will 
be made later.) Thus the j(x)'s are c-number sources and since Wick's theorem is used 
in perturbation theory, the field operators ip(x) are in the interaction picture. The normal 
ordered products iV[»] and contraction D[if)(x)ip(y)} have been defined above in terms of 
the arbitrary split of the fields into positive and negative parts; ip(x) = ip + (x) + tp~(x). 
Following the discussion in the previous section, we therefore make no assumptions about 
their symmetry properties. The first part of the proof presented here follows closely that of 
[3,4]. 

Recall that for two operators A and B 



if and only if 



e A+B = e A e B e ±[B,A] 



[A, [A,B]] = [B,[A,B}] = 0. 



(2.2) 
(2.3) 



Recall, too, that time ordering puts the larger times to the left, and so for example in the 
case of two operators A(ti) and A(t 2 ) with t 2 > ti, 



T 



^A(ti)+A(t 2 ) 



e A(t 2 ) e A( tl ) 



T 



e A( tl ) e A(t 2 ) 



. (twooporder) 



(2.4) 



We comment that a similar relation holds for any operation which places the two operators 
in a given order. This will be used for normal ordered products below in equation (2.11). 
Consider the LHS of (2.1), 



T 



3 — i J dt J d 3 x j(x)ip(x) 



T 



7 -i J dt o(t) 



(2.5) 



Let ti and tf denote the initial and final time. Divide the total time interval into iV equal 
intervals with At = t - 1 T p- and 



tf — tjy > ijv-i > . . . > ti > to — U 



By definition, 



[ tf dtO(t) 

hi 



N 



(2.6) 



(2.7) 



So here, 



T 


-i f tf dt 0(t) 
e Jt i 


= lim T 

N—>oo 










= lim 


e 


-iAt 0(t N ) e ~iAt 0(tjv_i) e ~ iAt °(*i) 



lim 



-iAf X£ =1 0(t k ) e -i(At) 2 E 1<k< i <N lO(tl),0(tk)] 



(2.8) 
(2.9) 
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In going from (2.8) to (2.9) we have used (2.2), which is legitimate since [0(ti),0(tk)] is a 
c-number. Note also that in this exponential t\ > t k . Thus using the definition of 0(t), (2.9) 
gives 



T 



7 -i J d 4 xj(x) i/>(x)l _ -i J d 4 xj(x) ip(x) -± J d 4 x J d 4 y 9(x -yo)j(x) [i/>(x) ,ij>(y)] j(y) 



(2.10) 



In order to relate the time ordered product in (2.10) to a normal ordered product, consider 
an arbitrary split of the fields into 'positive' and 'negative' parts introduced in (1.6) and 
the corresponding definition normal ordering given in the previous section. Using (2.2), the 
first exponential on the RHS of (2.10) may be re- written as 



-i J d 4 x j(x) tp(x) 

= e -i J d 4 xj(x) {i,-(x)+f+(x)} 

— g-» / d 4 xj(x) ip~(x) -i J d 4 xj(x) tp+(x) 

± 2 jd 4 x J d 4 yj(x) [i,-{x)^+{y)\ j( 



N 



-i J d 4 xj(x) tf(x)l \ J d 4 x J d 4 yj(x) [i,-(x),f+(y)} j(y) 



(2.11) 



In obtaining (2.11) we have used (2.4) for normal ordered products, and (2.2), which requires 
that 



= 0. 



(2.12) 



This condition is satisfied if ip + and ip~ are linear in annihilation and creation operators, 
and in turn it ensures that the contraction is a c-number. Substituting (2.11) into (2.10) 
gives 



T 



x e 



-i J d 4 x j(x) 4>{x) 



N 



-i J d 4 x j(x) ip(x) 



\jd 4 x J d 4 yj(x) {d(x () -yoM(x)My)] -[V>" {x)^+{y)\ }j 



(2.13) 



Our proof now diverges from that given in the usual texts [3,4] since we are considering the 
possibility of having an arbitrary (unsymmetric) split of the fields rather than a (symmetric) 
one into annihilation and creation operators. Look at the part in curly brackets in the second 
exponential of (2.13); it is 

0(^0 -Vo) [^(x),ip(y)} - [^(x),^ + (y) 

= ^ + (y),^-(x)]+6(x -y ){[^ + (x),^ + (y)] + [^(x),^ + (y) 



+ 



0(yo-x o ) [ip + (y),ip (a;)] - B(x - y Q ) {[^ + (y) , ip + (x) 



+ 



= D [if>(y)if>(x)] , 
where we have used the definition of the contraction in (1.10). 



(2.14) 
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Thus, now substitute (2.14) into the final exponential of (2.13), relabel x <-> y and use 
the fact that the j(x),j(y) are c-numbers to switch their positions. This enables (2.13) to 
be re-written as 



T 



x 



exp{— i J d A x j(x)tp(x)} — N exp{— i J d 4 x j(x)tp(x)} 



d 4 x d 4 y j(x)D[ip(x)tp(y)]j(y)}, 



(2.15) 



thus proving (2.1). This equation is also true for fermionic operators. There the sources are 
Grassmann variables, and the contraction, defined in (1.10), contains anti-commutators. 



III. WICK'S THEOREM FOR NON-SYMMETRIC PRODUCTS 

Wick's theorem is obtained by functionally differentiating (2.15) with respect to the 
sources j(x). Observe that we need only work with an even number of fields as the pertur- 
bation is about the expectation value of the fields, and hence the expectation value of an 
odd number of fields will always vanish. 

Functionally differentiate (2.15) with respect to the sources 2n times and then set every 
j(x a ) = j a to zero. Since time-ordered products of fields evaluated at the same time points 
are undefined, the indices a refer to different space-time points and are therefore all different. 
The LHS of (2.15) gives 



T 



Sji ■ ■ ■ 5j 2n 



exp{-ijip} 



= T 



ja=0 



-ijip) 



2n 



Sji---5j 2n (2n)! 
T 



( ai > a2 >---' a2n 

ix (2n)! 

perm {a} v ' 

= 7l j2 , 3 ,...,2n- 

In (3.2) the sum takes ai, . . . a 2n through all permutations of 1, ... , 2n; 



(3.1) 

(3.2) 
(3.3) 



perm {a} = 



12 3 

Qi CL2 O3 



2n 

0>2n 



(3.4) 



and p is the number of times one has to interchange pairs of fermions in moving from 
1, . . . , 2n to ai, . . . ci2n order. The last identity (going from (3.2) to (3.3)) follows from the 
fact that by definition time ordered products are symmetric in their indices, and that there 
are (2n)! permutations of the numbers 1, 2, . . . In. This type of manipulation will not hold 
for normal ordered products and contractions since they are not, in general, symmetric in 
their indices. 

Carrying out a similar operation for the RHS of (2.15) gives a general form of Wick's 
theorem for non-symmetric products; 



7l,2,...,2n = £ £ (-If 
m=0 perm {a} 



N. 



01,02, ...,a2m 



1 



(2m)! (n — m)\ 




a2j-l,0,2j 



(3.5) 
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In the traditional split of the fields into annihilation and creation operators, the normal 
ordered products and contraction Df\ are symmetric as was discussed in section 1. 

As a result, the sum over permutations on the RHS of (3.5) can easily be evaluated and thus 
(3.5) rapidly reduces to the well known form of Wick's theorem. For two bosonic fields this 
is 

T^ 2 = N^+Dfl (3.6) 

and for four bosonic fields 

T m = < 3 ,4 + ^°K2 + + + ^>g) + 

4>i ( °3 } + + DflDfl + £>W^J + (3.7) 

Higher order relations follow the same pattern. Note that in the above expressions we have 
chosen to arrange the indices on any one given normal ordered product and contraction such 
that the smaller indices are to the left of the larger ones. In the case of symmetric products, 
this ordering is of course not important. However, it will also be used below where, starting 
from (3.5), we show that the above expressions (equations (3.6) and (3.7)) hold even for 
non- symmetric products N and D, as long as the contraction is a c-number. We again 
consider bosonic fields. The results below can be extended to fermionic fields by inserting a 
minus sign every time the positions of two fermionic fields are interchanged. The following 
relations and lemmas are useful for the proof. 

A. Lemmas and relations 

Lemma 1 Two point normal ordered products and contractions are related by 

N 2 ,i - N 1>2 = D 1>2 - D 2>1 . (3.8) 



Proof: By definition of the contraction (1.9) 

- 2;1 [i>t^l] + ^{[^t^t] + [^2^t] + K^l"]} 
= tpt,^l] + [^2,^1 

= A 2jl - N lt2 , 



(3.9) 



where 9i t2 = 6{t\ — 1 2 ) and where the last line follows directly from the definition of normal 
ordering (1.7). 

Lemma 2 Let X and Y be arbitrary operators. Then if the contraction is a c-number, 

N [XfafoY] = N [XfofaY] + (L> 2 ,i - D lt2 )N [XY] . (3.10) 
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Proof: By definition 



Therefore 



provided ipf,ip 



N 



N 



N [XfafoY] = N [Xtpf^Y] + N [Xi/if^Y 
+N [Xipi^Y] + N [xfcfoY 



N 
+N 
+N 



Xij+ij+Y 



+ N 



+ N 



X 



Y 



+ n x vr,^ 



Y 



N [XfofaY] + (£> 2)1 - D lt2 )N [XY] , 



and 



(3.11) 



(3.12) 



■0! , ip 2 are c-numbers. This implies that both if) + and ip are linear 



in the creation and annihilation operators, and hence that the contraction is a c-number. 



Relation 3 The number of ways of filling p indistinguishable bags each with exactly two 
marbles labelled 1, 2, . . . , 2q is precisely 



[(2q)(2q - l)][(2q - 2)(2q - 3)] . . . \{2q - 2p + 2){2q - 2p + (2 " ,! 



'p\2P (2q-2p)\2Pp\ 



(3.13) 



Consequence 3.1 Let labels Oj's take values from 1, 2, 3, . . . 2q. The number of distinct 
configurations of p contractions D ai a2 D a3 a 4 - ■ ■ D a2pl a2p with unordered indices is given by 
(3.13). 

Consequence 3.2 For a product of p contractions and one normal ordered product, with 
a total number of 2q indices; N aiA2 ^^ A2q _ 2v D a2q 

-2p+l ,<12<j-2p+2 ' ' ' Att,_i,o2,» the number of dis- 
tinct configurations with unordered indices is again given by (3.13). 

Using these, we now show that (3.6), (3.7) and the usual higher order version of Wick's 
theorem holds even for non-symmetric normal ordered products and contractions, as long 
as the contraction is a c-number. 



B. Proof for n = 1 



This is the easiest case. For n—\, (3.5) gives 



Ti, 2 = -(iVi, 2 + JV 2i i + £>i, 2 + D 2 ,i) 



(3.14) 



Substituting (3.8) into (3.14) gives 
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Ti, 2 = 7Vi j2 + D li2 



(3.15) 



showing that for n — 1 (3.6) holds even for non-symmetric normal ordered products and 
contractions. (This result is of course consistent with the definition of the contraction; from 
(1.9), D h2 := T lj2 - N h2 and similarly D 2 ,i := T 2jl - N 2>1 = T h2 - N 2 ,i. Adding these two 
equations gives (3.14).) 



C. Proof for n = 2 

For n = 2, (3.5) gives 



2 



^1,2,3,4 - H 



N 1 / 1 

1 v ai,a2,...,a2m 1 



n ^ 



perm {a} m=0 [ ( 2 ™) ! (2 - 77l) ! \^. =m+1 2 



«2j-l ,«2j 



perm {a} 

4! 



A \ (-^01,02,03,04 "I" 6A^ aii(l2 Z) a3)a4 + 3D aia2 D a3a4 ) (3.16) 

perm {a} 

= ^y'S'ai, 02, o 3 ,a 4 (3-17) 

perm {a} 

where in the first term of (3.16) m = 2, in the second m — 1, in the third m = 0, and 
£01,02,03,04 w iH be defined later. Provided we choose S^i ,02,03,04 correctly (see below), proof 
that 7^2,3,4 is given by the usual expression (i.e. (3.7)) even for non-symmetric N and D will 
then follow directly from (3.17) if we can show that S aua2ta3ta4 is symmetric under interchange 
of any two indices. 

To do that, consider initially the second term of (3.16); 2~Z P erm{a} 6N aiiCl2 D a3 ^ a4 . The first 
important point to notice is that because of the sum over permutations, the order in which 
the cij indices are written is irrelevant; for example 

-^01,02-^03,04 — -^01,03-^02,04- (3.18) 

perm {a} perm {a} 

Secondly, the factor of 6 can be understood in the following way. From consequence 3.2 we 
know that for a product of one normal ordered product and one contraction with a total of 
2q = 4 indices, there are 4!/2!2 = 6 distinct configurations with unordered indices. Thus 
we can re- write the second term of (3.16) as a sum over the 6 distinct configurations with 
unordered indices: 

6-^01,02-^03,04 [-^01,02 -^03,04 -^01,03-^02,04 -^01,04-^02,03 

perm {a} perm {a} 

+-£^02,03-^01,04 + D a2;(l4 N ai;a3 + D a3t a 4 N ai >a2 \ . (3.19) 

Note that since the indices on the contractions and normal ordered products are unordered, 
there was still some freedom in writing down the above expression. We have made the 
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particular choice to write the Oj indices on each individual iV and each individual D with 
the smallest value of i to the left. The reason for this will become clear later. 

Now carry out a similar procedure for the third term of (3.16) involving the product of 
2 contractions: the 3 terms arise from the 3 distinct configurations of two contractions with 
unordered indices (see concequence 3.1), and we then chose to write the indices on each 
individual contraction such that the smallest value of i is to the left. Thus S aua2ta3ja4 in 
(3.17) is defined by 

Sa± ,02,03,04 ' -^01,02,03,04 -^01,02-^03,04 ~l" -^01,03-^02,04 -^01,04-^02,03 
D a2j a 3 N aij a 4 4~ D a2j(l4 N aia3 + D a3a4 N aij(l2 

~\~D ai !<12 D aza4 + D aia3 D a2a4 + D ai ±D a2as . (3.20) 

Hence is given by the RHS of (3.7) where instead of the symmetric 

normal ordered products and contractions one uses the more general non-symmetric normal 
ordered products and contractions: this was the reason why (3.16) was decomposed as above. 

To show that S^i ,02,03,04 is symmetric under the interchange of any two indices, consider 
a particular order of the indices, say a, b, c, d and then the effect of a <-> b say. Now 

S a ,b,c,d = N a>btC;d + D a>b N Cjd + D Cjd (N a>b + D a fr) 

+ [D a , c N bid + D biC N a>d ] + [D a4 N byC + D b4 N a , c ] 

+ [D a , c D b4 + D a4 D b , c \. (3.21) 

The pairs of terms in square brackets are invariant under interchange a <-> b. Now use (3.8) 
and (3.10) to show that the first line of (3.21) is also invariant under a <-> b: 

N a ,b,c,d + D a>b N c>d + D c>d (N a>b + D a , b ) 
= N b ^ c>d + (Z\o - D a , b )N c , d 

+D a:b N Cjd + D Cjd (N b:a + D b , a ) 
= N b ^ c4 + D b:a N Cjd + D c , d (N b:a + D b:a ). (3.22) 

Thus Safocj = S bj a tC;d . One can prove more generally that S a>b>c ^ d is invariant under the 
interchange of any two adjacent indices and hence under any of the 4! permutations of the 
indices. Thus from (3.18) and (3.21) it follows that 

Tl,2,3,4 = Wl,2,3,4 + £1,2^3,4 + + A,4^2,3 + ^2,3^1,4 + 

£>2,4^1,3 + ^3,4^1,2 + ^1,2^3,4 + ^1,3^2,4 + ^1,4^2,3- (3.23) 

Hence, both at the 2- and 4- point level, the usual form of Wick's theorem holds even if the 
contraction and normal ordered product are not symmetric. 

D. Proof for arbitrary n 

At the nth order, the proof is a generalisation of that given above for n = 2. However, 
before proceeding, we use some of the expressions encountered in the previous subsection to 
motivate the introduction of two new pieces of notation. 
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The first enables us to write the terms appearing in square brackets on the RHS of (3.19) 
in a much more compact way. Recall that in this square bracket (which we now denote by 
Q), the i subscript on each normal ordered product and on each contraction is arranged in 
increasing order going from left to right. The simplest case is when a x = 1, a 2 = 2, a 3 = 3 
and 04 = 4, then Q is given by 

Q = D lt2 N 3A + D 1>3 N 2A + D 1A N 2>3 + D 2 , 3 N 1A + D 2A N lj3 + D 3A N 1>2 . (3.24) 

However, we can re-write this as 

Q= £ D bub2 N b3M (3.25) 

perm {b}6-B 

where the sum over permutations of {b} is restricted to include only those permutations 
which are an element of the set B defined by 

B = {perm {b} : b x < b 2 , b 3 < b 4 }. (3.26) 

(It is possible to check that (3.25) indeed reduces to (3.24).) If we now return to the case 
in which the a^s are arbitrary then Q is given by 

Q= E D abi , ab2 N ab3 , ab4 (3.27) 

perm {b}e_B 

where B is defined as in (3.26). 

We now introduce the second piece of notation whose aim is to enable one to re-write 
expressions such as (3.27) in a slightly simpler way (with fewer subscripts). Define 

N(ai, a 2 , . . . , a 2n ) '■= AT aij02) ... ja2ii , 

D( ai ,a 2 ):=D auaa . (3.28) 

As an example of the use of this notation, use (3.27) and (3.28) to re-write (3.19) as 

E 6D(ai,a 2 ) N(a 3 ,a 4 ) = E S D(a bl ,a b2 ) N(a bs ,a b4 ). (3.29) 

perm {a} perm {a} perm {b}£B 

We now continue with the proof of Wick's theorem to nth order. 
• Step 1. 

For a given value of m, the RHS of (3.5) can be re-written as a sum over distinct 
configurations of iV's and D's, with the indices on each individual iV and each individual D 
arranged in increasing order going from left to right (compare with (3.29)). To do that, use 
Lemma (4): 

Lemma 4 The following identity holds 
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(2nV 11 
(2m)!(n-m)!2— ^ , ^ ° 2 ' ' " ' ' ° 2m) . II ^K-i> a 2j ) 

\ S \ ) perm {a} j=m+l 

n 

= Yl X! N(a 6l ,a 62 ,...,a 62m ) II D ( a b 2j -i,ab 2j ) 

perm {a} perm {b}GB j=rn+l 



where 



(3.30) 



perm {6} 



1 2 3 
&i b 2 h 



2n 

b 2n 



(3.31) 



and the restriction on the sum of permutations of {b} is such that 1) it arranges the indices 
on each individual normal ordered product and on each individual contraction in increasing 
order going from left to right, and 2) it orders the contractions by size of the leading element 
so as to avoid repetition o/ -Di, 2-03,4 . . . as D^^Di^ 



B = {perm {b} : 



&i < b 2 < . . . < b 2m 
b 2m +i < b 2rn +2 
b 2m +3 < b 2m+ 4 



ha-i < b 



'2n 



; b 2m+ i < &2m+3 < • • • < &2n-l}- 



(3.32) 



Proof: From Relation (3), for a product of n — m contractions and one normal ordered 
product with a total of 2n indices, there are (2n)!/(2m)!(n — m)\2 n ~ m distinct configura- 
tions with unordered indices. These distinct configurations are represented by the sum over 
permutations of the b labels. On a given contraction or normal ordered product, the order 
of the indices is arbitrary, but above the choice has been made by the restriction. 

• Step 2. 

Using Lemma (4), re-write the RHS of (3.5) as 



7l, 2 ,...,2n 
1 



\ ' ' perm {a} 



J2 N (a bl ,a b2 ,...,a b2m ) JJ D{a h2j _^a b2j ) 

m=Operm{b}GB j=m+l 



1 



v ' perm {a} 

This is a generalisation of (3.20) to arbitrary n. 



(3.33) 



• Step 3. 

Use Lemmas (1) and (2) to show that 
Lemma 5 For any two permutations, {a} and {a'} 

S({a}) = S({a>}). 



(3.34) 
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Proof: Consider two permutations related by {a} — > {a'} with 



a i+i 



o-i+i — a[ 

dj = a'j \/j ^ i, i + 1. (3.35) 

The only non-zero contributions to A := S({a}) — S({a'}) arise when both the indices Oj 
and a i+ i are either on the same normal ordered product, or on the same contraction. The 
reason is that when these indices occur on different normal ordered products or contractions, 
the restriction on the sum of permutations of {b} orders the terms in exactly the same 
way in both the {a} and {a'} permutation thus making A vanish. (These terms are the 
generalisation of those in square brackets in (3.21).) 
Thus we may write 

S({a}) - S({a'}) 

n I 

= J2 2 [ N (a bl ,a b2 ,...,a b2 J - N(a' bi ,a' b2 ,...,a' b2 J 

m=0 perm {b}eB 
n 

x I] D{a b2j _ 1 ,a b2j ) (3.36) 

j =777+1 

n II 

+ J2 J2 N(a bl ,a b2 ,...,a b2m ) 

m=0 perm{b}e_B 

77 



X 



Dfa, Oi+i) - D(a[, a- +1 )] JJ D(a b2j _ 1 ,a b2j ). (3.37) 



j=m+l 



Here the prime on the sum in (3.36) means that the sum over permutations of {b} is further 
restricted by the condition that + 1) G {bi,b 2 , ■ ■ ■ ,&2?n}- (Thus it guarantees that the 
indices and a^+i only occur on the and not on the Z)'s.) In (3.37), the double prime 
means that this sum is further restricted to have {i,i + 1} = {&2fc-i,&2fc} for some k > m, 
and this is the reason for which the contribution with (J = k) is excluded in the product. 
(Hence the double prime on the sum ensures that the indices a, and a i+ i both occur on a 
single D.) 

Noting that when m = there are no contributions to (3.36) (since it involves no normal 
ordered products), and using Lemma (2), (3.36) may be re-written as 

77 / 

X! X] N(a bl ,a b2 , . . . , Oj, a i+ i, . . . , a b2m ) - N(a' bl ,a' b2 , . . . , a-, a- +1 , . . . , a' b2m ) 

7Tt=0 perm {b}gB 

77 

x n D ( a b 2] -n a b 2j ) 

J =777+1 

77 77—1 

= S 51 N(a Cl ,a C2 , . . . ,a C2m _ 2 ) [D(a i+U ai) - D(ai,a i+1 )] Y[ D(a Caj _ 1 , a C2j ) 

777=1 perm {c}S-B j=?rt 

71 — 1 71 — 1 

= S ^V(a Cl ,...,a C2 J [^(a i+ i,ai) -D(aj,a m )] JJ D(a C2j _ 1 ,a C2j ) (3.38) 

p=0 perm {c}£-B J=P+1 
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where p = m — 1, Cj = bj for j < i and Cj = bj +2 for j > i. 

The second term of equation (3.37) vanishes for m = n (as then it contains no contrac- 
tions) and thus it can be re-written as 

n II 

J2 J2 N (a bl ,a b2 ,...,a b2m ) 

m=0 perm {b}G-B 



X 



D(a h Oi+i) - D(a[, a' i+1 j\ D(a b2j _ 1 ,a b2j ) 



j=m-\-l 



n-l 



= J2 J2 N(a bl) ... ) a b2m )[D(a i) a i+1 )-D(a i+l) a i )} 

m=0perm {b}GB 
n 

x J] D(a b2j _ i: a b2j ). (3.39) 

j — m + l 

Recalling the definition of the c's, it is now possible to see that (3.38) and (3.39) just differ 
by a sign, and thus from (3.37) we have proved that A = ^({a}) — S({a'}) = 0. Hence A 
vanishes for any two permutations. 

• Step 4. 

Since 5*({a}) = S({a'}) for any two permutations {a} and {a'}, then from (3.33) one has 

7l, 2 ,...,2n = 7^T| E S d a }) 
^ ' ' perm {a} 

= S({a}) (3.40) 
(wickns) 

where ^({a}) is just given by the usual form of Wick's theorem with the ordered indices. 
This finishes the proof that Wick's theorem holds in its usual form for non-symmetric normal 
ordered products and contractions, provided the contraction is a c-number. 

IV. CONCLUSIONS 

In this paper we have worked with a general split of the field operators into two parts 
and the corresponding non-symmetric normal ordered products and contractions. We have 
shown that, somewhat surprisingly, Wick's theorem still holds in its usual form for such non- 
symmetric products as long as the contraction is a c-number. As a result, Wick's theorem - 
an important calculational tool in field theory - is more general than was originally thought. 

We now comment on another interesting idea inspired by the discussion of Wick's theorem 
in the original paper on thermal field theory by Matsubara [9]. There, the definition of the 
contraction given for bosonic relativistic phonon field is always symmetric for all splits, and 
consistent with another definition of normal ordering 7V( T )[«]; 

A^[«2...Vg =JV[TfoM> 2 ...^„]]. (4.1) 
(nomat) 
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Note that N^[»] is distinct from N[»], the definition of normal ordering which was given 
in section 1 and was used through out this paper: by definition is always symmetric 
in the sense of (1.4) because of the way it is built out of the time-ordered product (hence 
our designation with a (T) superscript ' 12 1). It therefore gives a corresponding symmetric 
contraction 

£ (T) [</<i, ^2] = T 12 - = 0(ti - t 2 )[Tj>t,fc]* + <r0(*2 - hWt, Vila- (4.2) 

It is clear that there are in fact many possible splits of the field which may be used 
to define normal ordering, and even then the definition of shows that one can think 
of other definitions of normal ordering different from the N given in section 1. Whichever 
definition is used, we have shown that even if the normal ordered product is not symmetric 
and the contraction is a c-number, Wick's theorem can be used so that calculations are not 
too cumbersome. 

In practice the only definite cases we know of have symmetric normal ordered products. 
Symmetry breaking problems expressed in terms of Bogoliubov transformations [6,7] and 
equilibrium thermal field theory, whether using path-ordered methods [1] or Thermo Field 
Dynamics [13], all use our definition of N but have symmetric products. It is easy to see why 
symmetric products are so common; it is because the most useful normal ordered product is 
one which has zero expectation value. Thus taking the expectation value of Wick's theorem 
implies that the c-number contraction must then be symmetric (since it is then equal to the 
expectation value of the two-point time-ordered function). In turn the operator version of 
Wick's theorem now shows that the normal ordered operator must itself be symmetric. 

However, there is at least one situation where we might envisage the appearance of non- 
symmetric products, and that is non-equilibrium field theory. In this case the enviroment is 
changing with time so we should expect the split, and hence the normal ordered products, 
also to change in time. It is then not at all obvious whether one can ever define a normal 
ordered products of fields at arbitrary times which have zero physical expectation values. 
Thus, reversing the argument of the previous paragraph, it is not at all clear that symmetric 
products will simplify the calculations, and therefore it is reasonable that non-symmetric 
products may prove to be a useful tool. 
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